J Wood Sci (2013) 59:375-382 
DOI 10.1007/sl0086-013-1342-6 


ORIGINAL ARTICLE 


Effect of annual rings on transmission of 100 GHz millimeter 
waves through wood 

Soichi Tanaka • Yuko Fujiwara • Yoshihisa Fujii • 

Shogo Okumura • Hiroyoshi Togo • 

Naoya Kukutsu • Shoji Mochizuki 


Received: 25 December 2012/Accepted: 1 April 2013/Published online: 16 April 2013 
© The Japan Wood Research Society 2013 


Abstract For the application of millimeter wave (MMW) 
technique to nondestructive evaluation of wood, the effect 
of annual rings on the behavior of a 100 GHz MMW 
transmitted through wood was examined. The complex 
amplitude was measured for 2 mm thick flat- (LT) and 
quarter-sawn (LR) specimens of several species with dif¬ 
ferent annual ring structures at 11 % moisture content, of 
which the density distribution was measured using X-ray 
radiography. For the LT specimens of all species and the 
LR specimens with small density fluctuation or with ear- 
lywood width smaller than the wavelength of the MMW 
(=3 mm), the amplitude and phase of the transmitted wave 
were similar to those of the wave without a specimen. For 
the LR specimen with large density fluctuation and with 
earlywood width close to or larger than the wavelength, the 
amplitude and phase were different from those of the wave 
without a specimen. All the measured complex amplitudes 
were well expressed using a diffraction model. It was 
concluded that the MMW is deformed by the density dis¬ 
tribution, and then its components with periods shorter than 
the wavelength by diffraction are attenuated. 
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Introduction 

Wood is widely used in industrial products such as furni¬ 
ture and structural members. The constant use of wood in 
industry, however, has been sometimes disturbed by the 
fluctuation of physical properties in wood because this 
fluctuation may cause defects in the wood products. For 
example, the abrupt change of wood density causes the 
structural defects such as raised and loosened grains. Thus, 
the physical properties should be evaluated nondestruc- 
tively to control the quality of the wood products [1]. 

A millimeter wave (MMW) technique, which deals with 
electromagnetic waves in the 30-300 GHz frequency 
range, has considerable promise as a tool for the nonde¬ 
structive evaluation (NDE) of the wood products, because 
it is inexpensive, non-contact, non-invasive, and compact. 
Recent studies have shown that this technique can obtain 
the transmission image of wood [2, 3] with higher spatial 
resolution than the conventional microwave technique in 
the 3-30 GHz frequency range [4-13]. However, the image 
presented no quantitative physical property of wood, 
because the effect of wood structure on the MMW trans¬ 
mission through wood was obscure. 

The transmitted waves are probably affected by micro- 
and macrostructures [14]. The micro structure can be 
regarded as homogeneous because wood cells, such as the 
tracheids, fibers, and ray cells, have cross-sectional 
dimensions much smaller than the wavelength of the 
MMWs, 1-10 mm. The macrostructure such as annual 
rings can be regarded as inhomogeneous because annual 
ring widths are close to or longer than the wavelength. 
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The effect of the microstructure on the MMW trans¬ 
mission has been explained by dielectric properties of 
wood [15]. On the other hand, little attention has been paid 
to the effect of the macrostructure such as annual rings. 
The understanding of this effect may improve the mea¬ 
surement accuracy of the physical properties of wood. 

The purpose of this study was to examine the effect of 
annual rings on the behavior of the MMWs transmitted 
through 2 mm thick flat- and quarter-sawn specimens of 
softwood and ring- and diffuse-porous hardwood species. 
The complex amplitude of the transmitted wave was 
measured using linear scanning and discussed based 
on a model in consideration of both the micro- and 
macrostructures. 


Theory 

Macro structure such as the annual rings is formed by the 
aggregation of microstructure such as the wood cells. Thus, 
in this study, the effect of the macrostructure on the 
MMWs is treated as the distribution of the effect of the 
micro structure on them. 

When an MMW is transmitted through a wood specimen, 
it will be locally affected by the complex permittivity due to 
the microstructure and be diffracted in the specimen simul¬ 
taneously. However, it is difficult to build a model for such a 
phenomenon. Let us consider that the diffraction starts at the 
plane which is perpendicular to the wave propagation 
direction and includes the center of the specimen in the 
direction, while the wave is affected by the local complex 
permittivity in all the position in the specimen. Thus, it can 
be assumed in a right-handed Cartesian coordinate system 
(x, y, z) that a wave is rectilinearly propagated through a 
specimen in the z direction, while being locally affected by 
its complex permittivity and is subsequently diffracted in 
free space at a distance of half the optical length of the 
specimen along the z-axis in the specimen and at a distance of 
/ 0 (=1 mm) in air (Fig. 1) if the double counting of the phase 
was corrected. The wave was assumed to be polarized con¬ 
stantly and linearly in the y direction, to illuminate perpen¬ 
dicularly the specimen whose longitudinal direction is in 
the y direction (Fig. 1) and to be expressed as a complex 
amplitude. The complex permittivity of the specimen in the 
longitudinal direction was assumed to be constant along the 
z-axis and represented by s(x, y). 

The complex amplitude of the wave propagated rectilin¬ 
early through the specimen is formulated as follows [16,17]: 

f(x, y)=s 2 i (x,y)Vi(x,y), (1) 

where V\(x, y) represents the complex amplitude of the 
incident wave and S 2 i(x, y), the ^-parameter which is 
formulated as follows: 


Vi(x,y) f(x,y) 


g(x,y) 



propagation j n specimen in air 
in specimen 


Fig. 1 A diffraction model of MMW during propagation through a 
specimen. Vfa, y), f(x, y), and g(x, y), complex amplitudes of an 
incident, a rectilinearly propagated, and a diffracted waves, respec¬ 
tively; dRey/t/2 and Z 0 , distances of the diffraction areas in the 
specimen and air, respectively; e(x, y), complex permittivity distri¬ 
bution in the specimen; d and e, the thickness and averaged complex 
permittivity of the specimen, respectively 


l-{R(x,y)} 2 

S 2 i(x,y) = - -—- — - -A(x, y) 

1 - R{x,y)A{x,y) 


R(,, y) = ' ~ vpi 

1 + y/e(x,y) 


\ 


V 


A(x,j) = exp|i2Tt(y s(x,y) - Vjd/X^ j 



where d represents the specimen thickness (=2 mm) and 2 
the wavelength of the MMW in air (=3 mm). 

The diffraction of the waves in free space is explained 
by Fourier optics [18], in which the complex amplitude of a 
diffracted wave at output xy-plane is related to that of a 
wave at input xy-plane using a transfer function. When the 
complex amplitude of the wave at input xy-plane is equal to 
/(x, y) (Fig. 1), that of the diffracted wave, g(x, y), can be 
formulated as follows: 


oo oo 


g(x,y) = exp(-i2jt//A) J J H(v x ,v y )F(v x ,v y ) 


—oo —OO 


x exp{i27i(v x x -j- v >; y)}dv x dv } ;, 



where v x and v y represent the spatial frequencies in the x 
and y directions, respectively, exp(—i27i//2) represents the 
factor which prevents the double counting of the phase, 
and F(v x , v y ), H(v x , v y ), and / represent the Fourier 
transformation, transfer function of /(x, y), and the optical 
length between input and output planes, respectively, 
which are formulated as follows: 
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oo oo 


F{V X ,Vy) 


f(x,y) exp{—i 2n{v x x + v y y)}dxdy, 


—oo —oo 


H(v x , v y ) = exp<j i 2% { 2 _ v 2 _ v 2 ) / 


/ — /o + 


JRe\/£ 


(4) 


(5) 


( 6 ) 


where s represents the average complex permittivity of the 
specimen, 


s = 


f- 0 c f-oo <x,y)dxdy 
f 00 f°° cLvdv 

J — oo J — OO 


(7) 


Equations (1)—(7) define a diffraction model in this 
study and are summarized as follows: 


X-ray 
source 
r 32 kV 
. 1.0 mA 


Wood 
specimen 


n 



X 


Image 

intensifier (I. I.) 



Intensity of X-rays, I{X, F) 

(640 x 480 pixel, 

256 gray levels) 
for LR specimen of sugi 

Fig. 2 Experimental set-up for X-ray radiography of wood 
specimens 


1. The incident wave is deformed by the complex 
permittivity distribution (Eqs. 1,2). 

2. The diffraction attenuates the deformed wave compo¬ 
nents with periods shorter than the wavelength, 2 
(Eqs. 3-7), because the transfer function H(v x , v^) acts 
as an attenuator if the spatial period of f(x, y), 

( v x + v y) _1/2 , * s l ess ^ an 


Experimental set-up 


specimen were converted into a visible image by an image 
intensifier (I.I.). The image was captured by a CCD camera 
(640 x 480 pixel, 256 gray levels) and sent to a computer 
to obtain the intensity of the rays, /(X, F), in a Cartesian 
coordinate system (X, F), where X and F are in the width 
and height directions of the image, respectively (Fig. 2). 
The specimen was set so that its center was at (X, F) = 
(320, 240) and its longitudinal direction was parallel to the 
F-axis. 

The intensity of the rays that passed through the speci¬ 
men at (X, Y) is formulated as follows: 


Specimens 

Flat- (large surface parallel to LT plane, LT) and quarter- 
sawn (LR) specimens, 100 mm square and 2 mm thick, 
were prepared from sound wood of softwoods, hinoki 
(Chamaecyparis obtusa ), sugi (Cryptomeria japonica ), and 
akamatsu (Pinus densiflora ); ring-porous hardwoods, kiri 
(Paulownia tomentosa) and keyaki (.Zelkova serrata ); and 
diffuse-porous hardwood, tochinoki (Aescuius turbinata). 
The specimens were stored over a saturated NaBr solution 
in desiccators at 27.5 °C until their moisture content 
reached 11 %. A Cartesian coordinate system (x, y, z) was 
defined on the specimen with its larger surface parallel to 
the xy-plane and its center on the z-axis. 

Density measurement with X-ray radiography 

The density of specimens was measured using X-ray 
radiography. Figure 2 shows an experimental set-up for 
X-ray radiography using an X-ray television inspection 
equipment (FI-30, Shimadzu). The X-rays were generated 
by an X-ray tube at a tube voltage and current of 32 kV and 
1.0 mA, respectively. The rays that passed through the 


/(X, F) = 7 0 (X, F) exp {-pp(X, Y)d }, (8) 

where 7 0 (X, Y) represents intensity of the rays at (X, F) 
without a specimen, p(X, Y) the air-dry density of the 
specimen averaged along the z-axis, and p the mass 
attenuation coefficient. Assuming that p is constant for all 
the specimens, Eq. (8) can be transformed into: 

r = -npj, (9) 

where V and p a represent the averages of the logarithmic 
intensity attenuation and air-dry density over the coverage 
area of the image, respectively, and they are expressed 
using discrete values for each pixel as follows: 

^640 ^480 , I{XJ) 
p _ 2^X=l Z^Y=1 m l 0 (X,Y) 

640 x 480 

0 _ Et°! Egl p(x, Y ) 

The relation of the average attenuation, T, to the product 
of the average air-dry density and the thickness, p a d , was 
obtained for all specimens in a preliminary experiment 
(Fig. 3). There was a good linear relation between them, 
and the mass attenuation, p , was determined to be 


( 10 ) 

( 11 ) 
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0.4640 cnT/g using Eq. (9) from the slope of the regression 
line. 

The density of a specimen, p(X, F), was estimated using 
Eq. (8) for 640 x 480 pixel points and resampled to 
p(x, y ) with 1024 x 1024 pixels over an area of 
—24 mm < v < 24 mm and —20 mm < y < 20 mm by 
Shepard’s method [19] using the X-ray image of a square 
mesh of 0.35 mm thick steel wires with 5 mm spacing to 
perform a fast Fourier transform (FFT) and an inverse FFT 
(IFFT) in the discussion described below. Furthermore, an 
oven-dry density, p 0 (x, y), was obtained from the air-dry 
one, p(x, y), using the equation, p = 1.11 p 0 /( 1 + 0.094p o ), 
for wood at 11 % MC [20] . 

Oven-dry density, p 0 (x , y), obtained in this way is shown 
in Figs. 4 and 5 as grayscale images for the FT and FR 



Fig. 3 Relation of average attenuation, F, to product of average air- 
dry density, p a , and thickness, d, for all specimens 


specimens, respectively. The FT specimens showed 
smooth tones of gray (Fig. 4), while the FR specimens 
showed remarkable stripped patterns due to the annual 
rings (Fig. 5). 

The complex permittivity in the longitudinal direction of 
wood at 11 % MC, s(x, y), was determined as follows: 


s(x,y) = 1 + 


0.760(sws ~ l)p 0 (x,y) 
1 +0.094 p 0 (x,y) 



This equation was derived from a parallel capacitor 
formula, s = 1 — v ws + v ws &ws [21, 22], where s ws 
represents complex permittivity of wood substance and 
v ws , volume fraction of wood substance which was 
expressed as 0.760p o /(l + 0.094p o ) at 11 % MC [23]. In 
an experiment described in another article submitted to this 
journal, Eq. (12) was confirmed to be well fitted to the 
complex permittivity in the longitudinal direction for FT 
specimens of nine wood species at 11 % MC measured 
using a free space method in the range of p 0 = 0.2-1.0 g/ 
cm , and the complex value of e ws was estimated to be 4.1 
(1 + 0.12i). 


Apparatus and measurement 

Figure 6 shows an experimental set-up of an MMW mea¬ 
surement system arranged in a right-handed Cartesian 
coordinate system (v, y, z). An MMW signal of 100 GHz 
with an average output power of 10 mW was oscillated in a 



(c) Akamalsu 



(d) Kiri (e) Keyaki (f ) Tocfimoki 
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Fig. 4 Images of oven-dry density distributions in xy-plane for LT specimens 
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Fig. 5 Images of oven-dry density distributions in xy-plane for LR specimens 



Fig. 6 Experimental set-up of an MMW measurement system 
arranged in a right-handed Cartesian coordinate system (x, y, z). 1 
Gunn diode (100 GHz), 2 Isolator, 3 Directional coupler, 4 PIN diode, 
5 Phase shifter, 6 Conical horn antenna, 7. Teflon piano-parabolic 
lens, 8 Specimen, 9 Open-ended rectangular waveguide (inner cross- 
section: x x y = 2.54 x 1.27 mm) covered with millimeter wave 
absorbent except open end, 10 Quadrature phase shifter, 11 Switching 
frequency generator (100 kHz), 12 Mixer, 13 Lock-in amplifier, 14 
IQ-mixer, and 15 Dielectric waveguide. V { and Vq, output voltages of 
in-phase (I) and quadrature (Q) components of an MMW signal 


Gunn diode and divided into two signals, measurement 
and reference, by a directional coupler. The reference 
signal was sent to an in-phase/quadrature mixer 
(IQ-mixer) through a dielectric waveguide. The mea¬ 
surement signal was modulated using a PIN diode at a 
switching frequency of 100 kHz and passed through a 
phase shifter. The wave linearly polarized in the 
y direction was transmitted by a conical horn antenna 
toward the z direction and collimated using a Teflon 
piano-parabolic lens with a diameter of 80 mm and a 
focal length of 33 mm. The wave that passed through the 
specimen was received by an open-ended rectangular 
waveguide with an inner cross-section of 2.54 mm in the 
x direction and 1.27 mm in the y direction and sent to 
the IQ-mixer, where the in-phase (I) and quadrature 
(Q) components were discriminated by comparing the 
measurement and reference signals. The signals from the 
IQ-mixer were processed into the voltages, Vi and Vq, in 
a lock-in amplifier. The front face of the waveguide, 
excluding the aperture, was covered with MMW absor¬ 
bent material to prevent reflections from the waveguide 
and surrounding apparatuses. 

The complex amplitude of the MMW was determined 
as: 

_ . / Vi Vq 

V Vl +1 \tan(^Q - &) (A Q /Ai)sin(^ Q - &) 
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where A h Aq, and are the amplitude and phase of the 
sinusoidal functions obtained as the regression lines of the 
relation between V h V Q and the phase controlled by 
the phase shifter and were determined in a preliminary 
experiment. 

The waveguide was set at (y, z) = (0, 113 mm) and 
moved along the v-axis in steps of 1 mm within a range 
of —20 mm < x < 20 mm for line scanning with and 
without a specimen (Fig. 6). Furthermore, it was set at 
z = 110 mm for area scanning of an incident wave 
(—24 mm < x < 24 mm and —20 mm < y < 20 mm, 
1 mm step along the x- and y-axes) (Fig. 6). 

Figure 7a and b show the distributions of the amplitude, 
Wi(x, y)l, and phase, argVi(v, y), of the incident MMW at 
z = 110 mm. Figure 7a shows a bright region in the center 
and concentric bright rings around it, so called main 
lobe and side lobes, respectively [24-27]. In Fig. 7b, 
there are three concentric regions, gray, dark gray, 
and black, corresponding to one main and two side lobes. 
These distributions were obtained over the area of 
—24 mm < x < 24 mm and —20 mm < y < 20 mm with 
49 x 41 pixels. To perform the FFT and IFFT in the dis¬ 
cussion described below, they were resampled to the dis¬ 
tributions with 1024 x 1024 pixels using Shepard’s 
method [19] and the complex amplitude, V\(x, y), was 
obtained. 

Figure 7c shows the distributions of the amplitude, I VqI, 
and phase, argV 0 , of the MMW without a specimen along 
v-axis at (y, z) = (0, 113 mm). The distribution, IV 0 I, has 
a large peak in center and two small peaks near either 
end, and argV 0 has three plateau regions in the range 
of —12 mm < v < 12 mm,—20 mm < x < —15 mm, and 
15 mm < v < 20 mm, corresponding to one main and two 
side lobes. 

Results and discussion 

The amplitude, IVI, and phase, argV, of the transmitted 
MMWs measured along v-axis at (y, z) = (0, 113 mm) are 
shown in Figs. 8 and 9 for the LT and LR specimens, 
respectively. For the LT specimens of all species (Fig. 8) 
and LR specimens of hinoki, keyaki, and tochinoki 
(Fig. 9a, e, f), the variation patterns of IVI and argV were 
similar to those of IV 0 I and argV 0 in Fig. 7c, though IVI was 
slightly attenuated in comparison with IV 0 I and argV was 
shifted away from argV 0 toward the positive direction. For 
the LR specimens of sugi, akamatsu, and kiri (Fig. 9b-d), 
on the other hand, the patterns of IVI and argV were dif¬ 
ferent from those in IV 0 I and argV 0 (Fig. 7c). 

The complex amplitude of a diffraction model, g(x, y), 
was estimated for z = 113 mm by substituting the inter¬ 
polated complexes, Vi(v, y) and po(x, y), for Eqs. (l)-(7) 



-20 -10 0 10 20 
x [mm] 



Fig. 7 Images of a amplitude, IVJ, and b phase, argV i? of incident 
wave at z= 110 mm, and c distributions of amplitude, I VqI, and 
phase, argV 0 , of wave without a specimen along x-axis at (y, z) = 
(0, 113 mm) 

and (12), of which Eqs. (3), (4), and (7) were computed 
using the IFFT, FFT, and numerical integration, respec¬ 
tively, in a discrete-finite form. In consideration of the 
aperture area of the receiving waveguide in Fig. 6, the 
simple moving averaged complex amplitude, g s (v, y), was 
calculated as simple moving averages of g(x, y) over 
2.54 mm in the v direction and 1.27 mm in the y direction. 
The amplitude, lg s (x)l, and phase, argg s (v), estimated from 
the diffraction model are shown as solid and dotted lines, 
respectively, in Figs. 8 and 9. The plots, IVI and argV, for 
all the specimens matched well the lines estimated from the 
model (Figs. 8, 9). This indicates that the diffraction model 
is valid for MMWs transmitted through wood specimen 
and that the waves may be deformed by the complex per¬ 
mittivity distribution due to the density distribution as 
shown in Figs. 4 and 5, and their components with periods 
shorter than the wavelength of the MMWs, 2 (=3 mm), 
must be attenuated. 
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Fig. 8 Distributions of amplitude and phase of transmitted MMWs for LT specimens measured along x-axis at (y, z) = (0, 113 mm). V, 
measured complex amplitude; g s (x), complex amplitude estimated from a diffraction model 
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Fig. 9 Distributions of amplitude and phase of transmitted MMWs for LR specimens measured along x-axis at (y, z) = (0, 113 mm). V, 
measured complex amplitude; g s (x), complex amplitude estimated from a diffraction model 


The LT specimens of all species and LR specimens of 
hinoki, keyaki, and tochinoki had small density fluctuations 
in the xy-plane or early wood widths smaller than A (Figs. 4, 
5a, e, f), while the LR specimens of sugi, akamatsu, and 
kiri had large density fluctuations in the xy-plane and 


early wood widths close to or larger than A (Fig. 5b-d). This 
is why the amplitude and phase of the transmitted wave for 
LR specimens of sugi, akamatsu, and kiri were different 
from those of the wave without a specimen on the basis of 
the diffraction model. 
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Conclusion 

The complex amplitude was measured for 2 mm thick flat- 
(LT) and quarter-sawn (LR) wood specimens with different 
annual ring structures, of which the density distribution 
was measured using X-ray radiography. For the LT spec¬ 
imens of all species and the LR specimens with small 
density fluctuation or with earlywood width smaller than 
the wavelength of the MMWs, 2 (=3 mm), the amplitude 
and phase of the transmitted wave were similar to those of 
the wave without a specimen. For the LR specimen with 
large density fluctuation along its large surface and with 
earlywood width close to or larger than 2, the amplitude 
and phase were different from those of the wave without a 
specimen. All the measured complex amplitudes were well 
expressed using a diffraction model. It was concluded from 
these findings that the incident wave is deformed by the 
density distribution and then its components with periods 
shorter than 2 by the diffraction are attenuated. Further¬ 
more, these findings will not only improve the measure¬ 
ment accuracy of the physical properties of wood, but also 
provide a new technique such as the detection of the annual 
rings without scanning. Further study is necessary to 
examine these effects for the specimens whose large sur¬ 
faces are not parallel to LT and LR planes. 
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